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Fermions in nature come in several types: 
Dirac, Majorana and Weyl are theoretically 
thought to form a complete list. Even though 
Majorana and Weyl fermions have for decades re¬ 
mained experimentally elusive, condensed matter 
has recently emerged as fertile ground for their 
discovery as low energy excitations of realistic 
materials [1—6]. Here we show the existence of 
yet another particle - a new type of Weyl fermion 
- that emerges at the boundary between electron 
and hole pockets in a new type of Weyl semimetal 
phase of matter. This fermion was missed by 
Weyl in 1929 due to its breaking of the strin¬ 
gent Lorentz symmetry of high-energy physics. 
Lorentz invariance however is not present in con¬ 
densed matter physics, and we predict that an 
established material, WTe 2 , is an example of this 
novel type of topological semimetal hosting the 
new particle as a low energy excitation around 
a type-2 Weyl node. This node, although still a 
protected crossing, has an open, finite-density of 
states Fermi surface, likely resulting in a plethora 
physical properties very different from those of 
standard point-like Fermi surface Weyl points. 

Metals can exhibit non-trivial topological features [2]. 
Of these, the ones with vanishingly small density of states 
at the Fermi level, called semimetals (SM), stand out. 
For these, a distinction between topologically protected 
surface and bulk metallic states can still be made, and 
their Fermi surfaces (FS) exhibit a topological (T) char¬ 
acterization. Two kinds of TSM have attracted special 
attention: Dirac and Weyl SMs. In these materials a lin¬ 
ear crossing of two (Weyl) or four (Dirac) bands occurs 
at the Fermi level (see left panel of Fig. 1). The effective 
Hamiltonian for these crossings is given by the Weyl or 
gapless Dirac equation respectively. The Weyl crossings 
are protected from gapping, a statement of the massless 
nature of the Weyl fermion. In the following we limit the 
discussion to Weyl crossings only, although our results 
also hold for Dirac crossings. 

The appearance of Weyl points (WPs) is only possible 
if the product of parity and time-reversal (TR) is not a 
symmetry of the structure. When present, a WP acts as 
a topological charge - either a source or a sink of Berry 
curvature. A FS enclosing a WP has a well-defined Chern 
number, corresponding to this point’s topological charge. 
Since the net charge has to vanish in the entire Brillouin 
zone (BZ), WPs always come in pairs; WPs are stable to 



Figure 1. Possible types of Weyl semimetals. Left panel: 
Type I Weyl point with point-like Fermi surface. Right panel: 
Type II Weyl point is the touching point between electron and 
hole pockets. 

weak perturbations and can only be annihilated in pairs 
of opposite charge. A large number of unusual physical 
phenomena are associated with Weyl TSM, including the 
existence of open Fermi arcs on the surface FS [1, 7], and 
different magnetotransport anomalies [8-14]. 

Weyl SMs with broken TR symmetry have been pre¬ 
dicted to exist in several materials [1, 15, 16], but exper¬ 
imental verification has been lacking. More recently, the 
Weyl SM was predicted to exist in inversion-breaking sin¬ 
gle crystal nonmagnetic materials of the TaAs class [3, 4]; 
this prediction was soon verified experimentally [5, 6]. 

To present, Weyl SMs have been thought to have a 
point-like FS at the WP. We refer to these as type-I 
WPs (WPl), to contrast them with the new type-II WPs 
(WP2) that exist at the boundary of electron and hole 
pockets, as illustrated in the right panel of Fig. 1. We dis¬ 
cuss general conditions for WP2s to appear, and present 
evidence that WTe 2 , a material recently reported [17] to 
have the largest known to date never-saturating magne¬ 
toresistance, is an example of the the new type of TSM 
hosting eight WP2s. They come in two quartets located 
0.052 eV and 0.058 eV above the Fermi level. We present 
topological arguments proving the existence of the novel 
TSM phase in this material. Doping-driven topological 
Lifshits transitions, characteristic of WP2s, as well as 
emerging Fermi arcs at the surface, are presented. 

We start by considering the most general Hamiltonian 
describing a WP 
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where Aij is a 3x4 matrix, with i = and j = 

0,x,^,z being indices of the 2x2 unit and three Pauli 
matrices. The spectrum is given by 


e±(k) = ± 

i=l 


'i 


=T(k]±Uik) 


j = l i=l 


( 2 ) 

The kinetic part of energy, T(k), also linear in momen¬ 
tum, tilts the cone-like spectrum. Breaking the Lorentz 
invariance of the Dirac Weyl fermions in quantum field 
theory, this tilt was previously considered unimportant. 
Since Lorentz invariance does not need to be respected in 
condensed matter, its inclusion is important and leads to 
a finer classification of distinct FSs in a one-to-one corre¬ 
spondence with the theory of quadric surfaces. The con¬ 
clusion: there are exactly two distinct types of WPs [18] 
(see Supplementary Material). 

For a direction in reciprocal space, where T is domi¬ 
nant over t/, the tilt becomes large enough to produce 
touching electron and hole pockets, rather than a point¬ 
like FS. Thus, the condition for a WP to be of type II is 
that there exists a direction k, for which T(k) > f/(k). 
If such a direction does not exist, the WP is of Type 
I. The clear qualitative distinction in FSs of the two 
types of WPs leads to drastic differences in the thermo¬ 
dynamics and response to magnetic fields [19]. In partic¬ 
ular, in sharp contrast to WPl, which exhibits a chiral 
anomaly [8] for any direction of the magnetic field, the 
chiral anomaly appears in a WP2 when the direction of 
the magnetic field is within a cone where T(k) dominates. 
If the field direction is outside of this cone the Landau 
level spectrum is gapped and has no chiral zero mode [19] 
(see Supplementary Material). 

On the lattice, the no-go theorem [20] guarantees that 
Weyl fermions appear in pairs with opposite Chern num¬ 
bers. Since the Chern number of a WP is not changed by 
the kinetic term, WPs of different type can be chiral/anti- 
chiral partners of each other. The number of WPs of cer¬ 
tain type could be odd, but the total number of WPs has 
to be even (e.g. there can be one WPl and one WP2). 

We now proceed to describe WTe 2 , a material we iden¬ 
tified to host the new WPs. The crystal structure of 
WTe 2 is orthorhombic with space group Pnm2i (Cj^). 
Its primitive unit cell contains four formula units. The 
atomic structure is layered, with single layers of W sep¬ 
arated from each other by Te bilayers and stacked along 
the z-axis (see Supplementary Material). The distance 
between adjacent W atoms is considerably smaller along 
the x-axis, creating strong anisotropy. The unit cell has 
two reflection symmetries: mirror in the ^z-plane, and 
a glide plane formed by a reflection in the xz-plane fol¬ 
lowed by a translation by (0.5, 0,0.5). Combined they 
form a non-symmorphic two-fold rotation (^ 2 , which is 
important in the following symmetry arguments. 

The result of band structure calculations (see Supple¬ 
mentary Material) without spin-orbit coupling (SOC) is 
shown in Fig. 2(a) along the T-X direction, where an in¬ 



Figure 2. Band structure of WTe 2 . Panel (a): Band structure 
of WTe 2 without SOC. A fraction of T-X segment is shown: 
the point E has coordinates (0.375, 0, 0). Band gap of ImeV 
is shown in the inset signaling a gapless point nearby. Panel 
(b): Band structure of WTe 2 with SOC. Panel (c): One of the 
four pairs of WPs. Their locations are designated in reduced 
coordinates. Zero energy is set to Ef in all panels. 


termediate point E = (0.375, 0,0) is introduced. Apart 
from electron-hole pockets, 16 WPs per spin are found 
in WTe 2 in the absence of SOC (not shown in Fig. 2(a)). 
Half of these points occur at points of low symmetry 
at kz 7 ^ 0. The other half of SOC-free Weyls appear 
in the kz = 0 plane, where the product of TR and C 2 
{C 2 T = C 2 * T) forms a little group. Generically degen¬ 
eracies on high symmetry planes are forbidden; however, 
due to the C 2 T symmetry two-fold degeneracies are lo¬ 
cally stable at points in the kz = ^ plane [18]. On the 
T-X line the spectrum is generally gapped with a band 
gap of ~ ImeV, separating valence and conduction bands; 
see Fig. 2(a). 

Accounting for spin, but without SOC, bands become 
doubly degenerate corresponding to opposite spin projec¬ 
tions. This doubles the topological charge of each Weyl, 
since by SU(2) symmetry WPs corresponding to oppo¬ 
site spins have identical topological charge. Infinitesimal 
SOC cannot gap them, giving a general recipe to search 
for Weyl semimetals [18]. 

In WTe 2 , however, SOC is not small. When turned 
on, it preserves electron-hole pockets, but significantly 
changes the structure of WPs. At intermediate SOC WPs 
can move, emerge or annihilate in pairs of opposite chi¬ 
rality. At full SOC all WPs dX kz ^ are annihilated. 
In the kz = ^ plane, double degeneracies at isolated k- 
points are still allowed by symmetry [18]. Eight such 
gapless points are found, formed by the topmost valence 
and lowest conduction bands at full SOC. A pair of such 
points is shown in Fig. 2(c). The other three pairs are re¬ 
lated to this one by reflections. Energetically both points 
are located only slightly (0.052 and 0.058 eV) above the 
Eermi energy. See Supplementary Material for details. 

Establishing degeneracies of bands (WP) through ab- 
intio is prone to finite size effects: one can imagine an 
unfortunate scenario where a degeneracy point actually 
turns out to have a minuscule gap upon further increase 
in computing resolution. To rigorously establish the de¬ 
generacies of WPs, we have performed a multitude of 
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error-free tests that involve computing topological in¬ 
dices. These are presented in the Supplementary Ma¬ 
terial. We have computed the topological charge (±1) of 
each Weyl point found through an extension of the Will- 
son loop eigenvalue method [21, 22] to type II Weyls. 
We have then computed the Z 2 topological indices on 
a multitude of planes (including both standard and non¬ 
standard geometries) in the Brillouin zone. In total, these 
tests not only proved the existence of Weyl nodes but 
also were more than sufficient to elucidated the structure 
of the Berry flux connection between Weyls and of the 
Fermi arcs on the surface of WTe 2 . The Fermi arc struc¬ 
ture is consistent to the one presented through ab-initio 
calculations below. 

To check the new nature of the WPs, we obtained their 
dispersion from first-principles calculations and fitted it 
to the theoretical one derived by symmetry analysis (Sup¬ 
plementary Material). Keeping terms linear in the 
momentum relative to the position of the WP, we get: 

5k = Akx + Bky ± e^k‘1 -h {ak^ + ckyY + {bkx + dky)‘^ 

( 3 ) 

The values of the parameters ^4, B, a, 6, c, d and e are 
given in the Supplementary Material. The kinetic part 
of the energy dominates along the line connecting the 
nearest WP (see Fig. 2(c) ). We thus conclude that WTe 2 
is a type II Weyl SM. 

We now discuss the FS topology and possible topolog¬ 
ical Lifshits transitions in this compound. The evolution 
of the FS obtained from first-principles calculations is 
shown in Fig. 3 for different values of Ep. Due to reflec¬ 
tion symmetries, only part of the k^ = 0 cut of the FS is 
shown. For = 0 eV the FS is formed of two pairs of 
electron and two pairs of hole pockets (8 pockets total), 
separated in momentum space. For each pair a larger 
surface completely encloses the smaller one, in agreement 
with experiments [23]. This is shown in Fig. 3(a), where 4 
halved (2 n and 2 p) pockets are shown. The other halves 
are obtained by glide reflection pxz] remaining kx > 0 4 
pockets are obtained by mirror rriyz. All Fermi surfaces 
have zero Chern numbers in this case. 

When E^' is raised, two additional electron pockets 
appear; the previously existing electron pockets persist. 
The hole pockets shrink quickly, two disappearing com¬ 
pletely. Each of the other two split into two disconnected 
pockets. As a result, there are 6 electron- and 4 hole- 
pockets in total (see Fig. 3(b) for illustration). For the 
Fermi level tuned to the first WP Ep = 0.052 eV, (cor¬ 
responding to the addition of ~ 0.064 electrons per unit 
cell), each of the two newly appeared n-pockets touches 
two p-pockets at the positions of the WPs, as illustrated 
in Fig. 3(c) for the part of kz = 0 plane. Further increase 
of Ep disconnects the electron and hole pockets again 
- see Fig. 3(d) for Ep = 0.055 eV - but with changed 
topology: electron pockets still have zero Chern num¬ 
bers, since they enclose two WPs of opposite charge, 
related by pxz- The hole pockets have Chern numbers 
±1. Topologies of the other p-pockets can be obtained 



Figure 3. Fermi surface at = 0. A part of the BZ is 
shown, (a). Ef = OeV: electron (blue and green solid) and 
hole (red and magenta dashed) pockets come in pairs, (b). 
The structure of electron and hole pockets at higher energies 
{Ef = 0.055eV shown). There are four hole (I shown) and 6 
electron (halves of the 3 of them shown) pockets. The zoom 
of the framed region is shown in (c)-(e) for different values 
of Ef. (c). Ef = 0.052eV is set to the lower energy Weyl 
point. A touching between n- and p- pockets occurs at the 
Weyl point, (d). Ef = 0.055eV is set to be in between the 
two Weyl points. The n- and p-pockets are disconnected. 
The hole pocket encloses a C = -hi Weyl point. The electron 
pocket encloses the C = — I point and its mirror image (not 
shown); the net Chern number of this pocket is zero. (e). At 
the higher energy Weyl point Ef = 0.058eV electron and hole 
pockets touch again (shown). They reopen at larger Ef with 
zero Chern numbers. 


by Chern-number flipping mirror and glide symmetries. 
The pockets touch again (see Fig. 3(e)) when the Fermi 
level is tuned to the higher energy WP (Ep = 0.058 eV, 
with 0.079 additional electrons per unit cell). Upon 
raising Ef further, the pockets disconnect again, and all 
FS Chern numbers become zero. 

WTe 2 is close to a topological transition occurring 
when the closeby WP of opposite chirality annihilate each 
other. To facilitate the observation of topological Lifshits 
transitions, hydrostatic pressure can be applied. Neigh¬ 
boring WPs are pushed away from each other in /c-space 
under compression. In particular, a 0.5% (2%) compres¬ 
sion increases the distance between the WPs from 0.7% 
to 2% ( 4%) of the reciprocal vector IG 2 1 (see Supplemen¬ 
tary Material for a thorough discussion of strain effects, 
including how to obtain only four Weyl points). 

Finally, we discuss the topological surface states of 
WTe 2 . Due to reflection symmetries, WPs of opposite 
chirality are projected on top of each other on the (100) 
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Figure 4. Topological surface states. Left panel: spectral 
function of the (001) surface. Fermi level is set in between 
the Weyl points (shown in green). Right panel: (001) surface 
Fermi surface and a Fermi arc connecting p- and n-pockets. 
Green crosses mark the positions of Weyl points. 


Lifschitz transition occurs and the system acquires finite 
density of states at the Weyl node. Although several pos¬ 
sibilities of small tilting of the Weyl cone were discussed 
in some works [25-29], none of these works pointed out 
the Lifschitz transition that gives rise to a fundamen¬ 
tally new type of Weyl fermion, with new thermodynamic 
properties. Given the existence of at least two types of 
Weyl fermions, we can ask if any other types of Weyl 
points can exist. In the present section we show that 
a mathematical theorem states that only two types of 
Weyl points are possible, and hence our classification is 
complete. 

The most general linear in k Hamiltonian of a Weyl 
point is 


and (010) surfaces which hence do not exhibit topologi¬ 
cally protected surface states. For the (001) surface all 
the WPs project onto distinct points, hence topological 
surface states appear. When E’p is tuned in between the 
WPs, the hole pocket has non-zero Chern number and a 
Fermi Arc emerges from it, connecting it to the WP of 
opposite Chern number inside the n-pocket. Fig. 4 illus¬ 
trates the (001) surface spectral function where surface 
states connecting electron and hole bands are clearly vis¬ 
ible in the left panel. The fermi surface of these surface 
states is a topological Fermi arc (right panel of Fig. 4) 
connecting the topological p-pocket (Fig. 3(c-e)) to the 
n-pocket. The other surface state crossing the p-pocket 
emerges from the n-pocket (not seen in the figure) and 
goes back into it, and thus can be pushed into the con¬ 
tinuum of bulk states not to be seen at this Ep (see Sup¬ 
plementary Material). 

Of other transition metal dichalcogenides another 
strong candidate material is MoTe 2 [24] - reported to be 
a semimetal resembling pressurized WTe 2 . These can be 
used to explore new physical phenomena arising in this 
novel TSM phase. 
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Appendix A: Classification of Weyl points 

Type II Weyl points occur when type I Weyl nodes 
are tilted enough - along some specific direction - that a 
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i7(k) = V • k + ^ kiAijdj (Al) 

bi=i 

where A is a 3x3 matrix and cr’s are the Pauli matrices. 
The Eermi surface for the chemical potential ji can be 
obtained by solving the equations e±(k) = /i, where 


ej-(k) = v-k± y^^ki[AA'^]ijkj (A2) 

V ij 

are the two eigenvalues of the Hamiltonian. By squaring 
the Eq. A2 the Eermi surface can be shown to be sat¬ 
isfied on the quadric surface described by the following 
equation 

ki{[AA^]ij — ViVj)kj -f 2/iv • k — /i^ = 0 (A3) 

Classification of all possible quadric surfaces is 
known [30]. Applying this classification to the equation 
above one can show that there only two possible types of 
Weyl points: type-I with a closed point-like Eermi sur¬ 
face, and type-H with an open Eermi surface, as described 
in this paper. The classification [30] also gives rise to dif¬ 
ferent kinds of nodal lines and surfaces, which will be the 
subject of further study [18]. 


Appendix B: Transport signatures of type-II Weyl 
semimetals 

Here we show that the transport properties of the two 
different types of Weyl points are very distinct. In par¬ 
ticular, the type-H Weyl points give rise to a new kind 
of chiral anomaly. Here we report the main conclusions, 
and a more detailed explanation will be presented else¬ 
where [19]. We also discuss another difference between 
the two types of Weyl points, which appears in the den¬ 
sity of states, leading to very different thermodynamic 
properties. 
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1. Novel chiral anomaly in type-II Weyl semimetals 

Let us consider the Hamiltonian 

H{k) = Ck,^k-a (Bl) 

which realizes a Weyl point with a +1 Chern number. 
The type of this Weyl point depends on the value of the 
parameter C, being of type-II for \C\ > 1 and of type-I 
otherwise. This Hamiltonian is a much simplified version 
for that of a type-II Weyl, but has the advantage of being 
analytically tractable in some limits and will be used for 
qualitative purposes. 

If a magnetic field of magnitude B is applied along the 
^-direction the spectrum of the Landau levels is given by 

= Ch ± (B2) 

where i = {eBlc)~i is the magnetic lengthscale, e is the 
charge of electron, c is the speed of light, and n > 0. The 
0-th Landau level has the energy 

Eo = {C^l)k, (B3) 

and is unpaired (chiral). The sign of the electron velocity 
associated with the 0-th Landau level is given by the sign 
of (C + 1), as explained below. 

If the electric field is applied in the same direction as 
the magnetic field, chiral electrons start flowing in the 
^-direction. Since Weyl points come in pairs of oppo¬ 
site Chern numbers, the low-energy spectrum of a Weyl 
semimetal will also host an antichiral 0-th Landau level at 
the position of the opposite Chern number Weyl point. 
The current then appears flowing from one Weyl point 
to another [2, 9-11, 31-37]. This is a realization of the 
Adler-Bell-Jackiw chiral anomaly [38, 39] of the quan¬ 
tum field theory in the lattice system [8] . Observation of 
this effect was reported recently in Dirac semimetals and 
in the TaAs type-I Weyl semimetal family [12-14]. For 
the type-II Weyl semimetal similar behavior is expected 
for the field applied in the direction of dominating kinetic 
term, when the spurious large trivial almost compensated 
Fermi surfaces that exist in the material are discounted 
and the Weyl contribution is singled out. 

A drastic difference between the two types of Weyl 
points described by Eq. Bl arises when the magnetic 
field is applied along the x-direction, that is in the di¬ 
rection, in which the kinetic term does not dominate 
in the case of |(7| > 1. For the type-I the above dis¬ 
cussed chiral anomaly still appears. For the type-II case 
of |C| > 1, however, the Landau level spectrum becomes 
gapped and no chiral anomaly can be observed, as illus¬ 
trated in Figs. 5-6 =0.1 in all the illustrations). 

This is easily understood in the limit C ^ oo by solv¬ 
ing the H = Ckz hamiltonian in a B\\x field. In a more 
general case, a long analytic calculation [19] proves the 
absence of a chiral mode for B\\x and C > 1. In fact. 


£ . — .. E^'- 


kz kx 


Figure 5. Landau level spectrum for the type-I Weyl point 
of Eq. Bl with C = 0.5. Left panel: the magnetic field is 
applied along the ^-direction. Right panel: the magnetic field 
is applied in the —x direction. The chiral 0-th Landau level 
is present in both cases. 


Figure 6. Landau level spectrum for the type-II Weyl point of 
Eq. Bl with C — 2. Left panel: the magnetic field is applied 
along the ^-direction. The chiral 0-th Landau level is present. 
Right panel: the magnetic field is applied in the —x direction. 
No chiral Landau levels exist. 


when the field direction is rotated from z to x the chiral 
Landau level persists until the angle 0 between the field 
the z-axis, becomes such that cos^ 6> = 1/C^. Eor the 
directions, where cos^ 6> < 1/C^ the chiral Landau level 
disappears. 

This situation is generic in type-II Weyl points. The 
chiral Landau level appears only when the magnetic field 
is applied in the direction, along which the kinetic term 
is dominant. If the field is applied in the other direc¬ 
tion - no chiral Landau level appears. This dependence 
of the presence or absence of the chiral anomaly on the 
direction of the applied magnetic field is the main dif¬ 
ference between the two types of Weyl semimetals, with 
immediate consequences for the transport properties. 

There is one more important feature of the type-II 
Weyl points that does not appear in the type-I case. This 
is the possibility for the velocity of the chiral Landau 
level to be locally in /c-space different from the chirality 
dictated by the Chern number of the Weyl node. As an 
example, consider the Weyl point of the Hamiltonian Bl. 
The Chern number of this point is positive, however, for 
C < —1 the velocity of the chiral Landau level appearing 
in the i-field becomes negative as illustrated in the left 
panel of Eig. 7. 

However, this anti-chiral effect is a feature of the lin- 
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Figure 7. Chirality change at small kz- Left panel: Landau 
level spectrum for the type-II Weyl point of Eq. B1 with C = 
—2. The 0th Landau level appears to be a left mover for a 
Weyl point with Chern number +1. Right panel: Schematic 
plot based on on a possible realization of the closure of the 
Fermi surface. The paradox is resolved, when higher-order 
in k terms are added to the Hamiltonian to close the Fermi 
surface. These higher order terms render the chirality of the 
full region to +1. The chiral band is a right-mover with the 
velocity changing twice at small values of kz- The framed 
region is that, where chirality is changing around small kz. 

earized Hamiltonian. In a real semimetal higher order 
terms must exist, responsible for closing the Fermi sur¬ 
face of the type-II Weyl point. These terms will correct 
the chirality of the Landau level far from the Weyl point 
as illustrated in the schematic plot in the right panel of 
Fig. 7. 



Figure 8. The density of states of the pair of Weyl points in 
WTe 2 . The contribution of the rest of the Fermi surface is 
not included. 



Figure 9. Left panel: crystal structure of WTe 2 . Right panel: 
2. Thermodynamics of type-II Weyl points Brillouin zone of the rhombohedral unit cell. 


Another significant distinction of the two types of Weyl 
semimetals appears in the thermodynamic properties. 
For the type-I Weyl point, described by the Hamiltonian 
of Eq. B1 the density of states g{E) behaves like 

aiE) oc 

For the type-II this behavior changes and the correspond¬ 
ing density of states becomes 

m « p (f>2 - (B5) 

The constant term po arises due to the presence of un¬ 
bounded (in the linearized model) electron-hole pockets, 
and it depends on the cutoff. Due to this term the den¬ 
sity of states always remains positive. These formulas 
were verified versus numerical simulation of the Hamil¬ 
tonian Bl. 

To illustrate this claim further we plot the density of 
states that arises due to the pair of Weyl points in WTe 2 
in Fig. 8. The contribution of the rest of the Fermi sur¬ 
face to the density of states is not included. The two 
clearly visible parabolic peaks in the density of states are 
due to the two type-II Weyl points is. The difference in 
the energy position of the peaks and the energy of Weyl 


points reported in the main text is due to the difference 
between the tight-binding model used here with the ac¬ 
tual first-principles calculation. The peaks correspond to 
the energies of Weyl points in the tight-binding model 
used to generate Fig. 8. 


Appendix C: Crystal structure and computation 
details 

The crystal structure of WTe 2 and its Brilloin zone 
are illustrated in Fig. 9. Several crystal structures are 
reported in the databases [40, 41]. We use the latest 
one from the work of Ref. [41] that was obtained at the 
lowest temperature. That work explicitly explains mi¬ 
nor differences with the earlier works [40] by the higher 
temperatures in those measurements. 

For completeness, we list the experimental structural 
parameters taken from Ref. [41] that were used in this 
work. The lattice constants are a = 3.477, b = 6.249 and 
c = 14.018, and both W and Te atoms occupy 2a Wyckoff 
positions corresponding to ( 0 ,^, 2 :) and (1/2, 2 : +1/2). 
Values of x and y for this structure are listed in Table 1. 

Band structure calculations (both with and without 
spin-orbit coupling (SOC)) were performed in VASP [42] 
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W(l) 

W(2) 

Te(l) 

Te(2) 

Te(3) 

Te(4) 

y 

0.60062 

0.03980 

0.85761 

0.64631 

0.29845 

0.20722 

z 

0.5 

0.01522 

0.65525 

0.11112 

0.85983 

0.40387 


Table L Positions of atoms in the unit cell of WTe2 given 
as coordinates for Wyckoff positions 2a of the Pnm2i space 
group. Bracketed numbers following the element symbol indi¬ 
cate distinct Wyckoff positions. There are two distinct Wyck¬ 
off positions for W and four for Te atoms. 



Figure 10. Band structure of WTe2 without SOC. Ef is set 
to 0. 


In addition, it can be shown that without SOC no band 
crossings can generally occur along the T — X axis. In¬ 
deed, ab initio results show that a small band gap exists 
in the band structure along T — X, as illustrated in the 
inset of Fig. 10. 

First-principles calculations support our conclusions. 
Without SOC we find 8 Weyl points (not accounting for 
spin) in the kz = 0 plane. In addition, we find 8 more 
spinless Weyls at low symmetry points ki 0 for all i. 
The coordinates of Weyl points and their Chern numbers 
are given in Tab. IT 



kx 

ky 

fc. 

C 

1 

0.1054 

0.0087 

0 

+1 

2 

0.1742 

0.1200 

0 

+1 

3 

0.14992 

0.07845 

0.2329936 

-1 

4 

0.14992 

0.07845 

-0.2329936 

-1 


Table II. Coordinates and Chern numbers of 4 out of 16 
Weyl points appearing without SOC. The other 12 points 
are related by the ones listed by reflections being located at 
{—kx,ky,kz), {kx, —ky,kz) with Chern numbers —C and at 
{—kx,—ky,kz) with Chern number C. Points at 7^ 0 are 
symmetric about the xy-plane due to the compound symme¬ 
try formed by time reversal and two-fold rotation. 


ab initio code using PAW [43, 44] pseudopotentials with 
6s^bd^ and valence electron configurations for W 

and Te correspondingly. The PBE [45] approximation 
was used. Spin-orbit coupling was implemented in pseu¬ 
dopotentials. The energy cutoff was taken to be 260eV. 
Gaussian smearing of width 0.05eV and a 12 x 10 x 6 
F-centered /c-point mesh were used to perform Brillouin 
zone integrations. The conclusions of the paper were also 
verified using a more elaborate pseudopotential for W 
with semicore p-states included in the valence, that is 
with bp^Gs^bd^ valence electron configuration. 


Appendix D: Electronic structure without spin-orbit 
coupling 

The band structure of WTe 2 in the absence of SOC is 
shown in Fig. 10. To understand the nature of crossings 
in this case we carried out an extensive symmetry anal¬ 
ysis that will be reported elsewhere [18]. Here we only 
state the main conclusions. Using a two band model 

Hikx-) ^y) ^z) ^y) ^z) 4 “ ^ ^ ^y) 

i=x,y,z 

(Dl) 

where ai are the Pauli matrices, one can show that on 
the kz = 0 plane only two of the three di coefficients 
are linearly independent. The degeneracy occurs if for 
some {kx^ky) both of these coefficients vanish. Thus, 
there are two independent constraints on two functions 
of kx^ky^ the codimension is zero, and the solution is in 
general possible at isolated points in the kz = 0 plane. 


Now, when spin is taken into consideration but SOC 
is still neglected, each Weyl point becomes doubled, with 
doubled chirality due to SU(2) symmetry. When SOC 
is gradually increased the Weyl points split and start 
moving in the Brillouin zone. For a material with weak 
SOC, one can expect this motion to be the only effect of 
SOC, so that the Weyl points do not vanish in general. 
In WTe 2 , however, SOC is large. To see its effect on 
the structure of Weyl points in the BZ, we carried out a 
smooth interpolation between no SOC and full SOC band 
structures. In the process, other Weyl points appear and 
disappear (by merging points with opposite chiralities), 
so that the final arrangement of Weyl points reported 
in the main text cannot be understood in terms of rear¬ 
ranging the Weyl points that are present in the absence 
of SOC. 


Appendix E: Electronic structure with spin-orbit 
coupling 

The fully relativistic band structure is illustrated in 
Fig. 11. Symmetry considerations indicate the possibility 
to have locally stable point-like degeneracies in the kz = 0 
plane, where a little group formed by a combination of a 
two-fold rotation and time-reversal exists [18]. As will be 
explained elsewhere in detail, this suggests that double 
degeneracies can occur at isolated points in the kz = 0 
plane. 

In accord with this argument, eight Weyl points (listed 
in Tab. Ill), described in the main text, are indeed lo- 
































Figure 11. Band structure of WTe2 with SOC taken into 
account. Ey is set to 0. 



kx 

ky 

kz 

c 

1 

0.12184 

0.03825 

0 

+1 

2 

0.12141 

0.0454 

0 

-1 


Table III. Coordinates and Chern numbers of 2 out of 8 Weyl 
points appearing with SOC. The other 6 points are related by 
the ones listed by reflections, being located at {—kx,ky,kz), 
{kx, —ky^kz) with Chern numbers —C. 


Gated in the kz = 0 plane. A general form of the Hamil¬ 
tonian around a Weyl point in the kz = 0 plane can be 
constructed using the following two dimensional repre¬ 
sentation for the combination of C 2 and time-reversal 


Cot — —icJzK 


(El) 


where K is complex conjugation. This form of the sym¬ 
metry can be shown to be consistent with the representa¬ 
tions of other symmetries in this non-symmorphic space 
group. 

Keeping only terms linear in k and subjecting the 
Hamiltonian to this symmetry we find the general form 
of the Hamiltonian around a Weyl point in kz = 0 with 
SOC 

iT(k) = €(k)-{-{akx-\-cky)ay-\-{bkx-\-dky)az-\-ekz(Jx (E2) 
where 


e(k) = Akx + Bky 


(E3) 


For the Weyl point to be of type H, the kinetic part of 
the energy should dominate over the potential one in at 
least some direction in /c-space. To find such directions 
for the points in question, we plotted the ratio 


R = 


{Akx + Bky)‘^ 


e^kl -h {akx + ckyY + {bkx + dkyY 


(E4) 


on the circles defined by k‘^^ky = 10“^, where kx and ky 
are taken in reduced coordinates, and the constants are 
explained in the main text. The results are illustrated 
in Fig. 12. The region, where > 1, meaning that the 



Figure 12. The ratio R of Eq. E4 for the circle of radius 10“^ 
drawn around the two Weyl points. Red color corresponds 
to R > 1, while blue is R < 1. The arrows show the di¬ 
rection from one of the two points to the other. Panel (a): 
Weyl point at (0.12184,0.03825,0). Panel (b): Weyl point at 
(0.12141,0.0454,0). 


kinetic energy dominates and the Weyl points are of type 
n, is shown in red. 

Einally, we fitted the Hamiltonian of Eq. E2 to the 
first-principles band structure around the Weyl points to 
get the values of its parameters. They are provided in 
Tab. IV. 


point 

A 

B 

a 

b 

c 

d 

e 

1 

-2.739 

0.612 

0.987 

1.107 

0.0 

0.270 

0.184 

2 

1.204 

0.686 

-1.159 

1.046 

0.0 

0.055 

0.237 


Table IV. Parameter values for the Hamiltonian of Eq. E2 
around the two Weyl ponts given in [eVA]. Point 1 refers 
to the Weyl at (0.12184,0.03825,0), while point 2 is 
(0.1241,0.0454,0). 


Appendix F: Weyl Nodes and Topological Indices 

To rigorously establish the degeneracies of WPs, we ex¬ 
amine the structure of the Berry curvature *F(k) of Bloch 
states around the gapless points. Since a WP represents 
a monopole in Berry curvature, the flux of R through a 
surface enclosing it in /c-space is quantized to the total 
topological charge enclosed, and the ES of WPl, being 
a closed surface, has non-zero Chern number. In case 
of WP2, however, the ES is open, and cannot be used 
to compute its topological charge. Instead, we integrate 
the Berry curvature of N bands, where N is the num¬ 
ber of electrons per unit cell. A surface, on which these 
N bands have an energy gap to the higher energy bands 
and which encloses the WP2 can always be found and 
is used to compute the topological charge of a WP2, as 
illustrated below. 

Eor WTe 2 these surfaces were taken to be spheres en¬ 
closing the Weyl points; see ii and £2 in Eig. 13(a). Eor 
the sphere £i the Chern number is found to be C = — 1, 
proving the existence of a WP of negative chirality inside 
it. The sphere £ 2 ^ enclosing the second WP, is found to 
have (7 = +1. This calculation combined with the use of 
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Figure 13. Weyl points of WTe2. Panel (a): schematic illus¬ 
tration of the BZ cross-section at = 0 and Weyl points in 
it, their chirality marked with red and blue color, correspond¬ 
ing to C = +1 and C — —1. The red line corresponds to the 
BZ segment shown in Fig. 2(c). Panel (b): Wannier charge 
centers (Wilson loop) on the plane, shown schematically as 
contour C in panel (a). It corresponds to a non-trivial Z2 
invariant. 




Figure 14. Panel (a): schematic illustration of the integration 
paths used to calculate topological charges of Weyl points. 
Panel (b): motion of the center of charge around the sphere 
is shown schematically for in panel (a) of Fig. 3 of the main 
text. The Chern number of the enclosed Weyl point is equal 
to —1. 


symmetries maps out the structure of the eight WPs in 
the BZ of WTe 2 . They are sketched in Fig. 13(a). 

To understand which pairs of Weyl nodes are con¬ 
nected by the Berry flux, thus elucidating the structure 
of the surface Fermi arcs, several topological invariants 
are computed from first-principles calculations. TR al¬ 
lows the computation of the Z 2 -invariant [46] on TR- 
symmetric planes in the BZ, defined by ki = 0,7r, where 
i = z. One such plane, the kz = 0 one, hosts WPs, 
and does not have a well-defined invariant. For the other 
five planes, where the N = 72 bands are gapped from the 
band N the topological invariant is computed using 
a hybrid Wannier centers technique [21] and found to be 
trivial. The corresponding Wilson loops [22] are gapped. 

A non-trivial TR Z 2 topological invariant can, how¬ 
ever, be defined for WTe 2 . It describes the motion of 
Wannier charge centers on a curved surface that crosses 
neighboring WPs, as shown by contour C in Fig. 13(a) 
where we plot its projection onto the {k^, /Cy)-plane - the 
surface extends through the BZ in the /C; 2 -direction. If 
the curve passes through a point {kx^ky) then it passes 
through its TR-image {—kx^—ky) as well and the surface 
satisfies the conditions of the Z 2 -pump [47]. We calculate 
Z 2 invariant on this surface directly from first-principles 
calculations using the Z2Pack software [48]. Fig. 13(d) 
shows that pumping in this plane occurs between two TR 
invariant momenta: kx = ky = —tt (S) and kx = ky = 0 
(F). Thus, the plane (C^kz) exhibits a quantum spin 
Hall effect. If an open surface in the ^-direction is intro¬ 
duced, helical surface states appear along the curve C in 
the surface BZ, becoming part of the Fermi arcs. As the 
kx,y = O 5 kz planes have gapped Wilson loops, we con¬ 
clude that the Berry flux connects pairs of neighboring 
WPs, as illustrated in Fig. 13(a). This calculation gives 
another proof that WPs exist in WTe 2 , resolving their 
connectivity and the location of Fermi arcs. 


1. Computation of chiralities of type-II Weyl points 

The Fermi surface of a type-II Weyl point is open, so 
it cannot be used for integrating the Berry curvature in 
chirality computation. Instead, the Berry curvature of 
N bands, where N is the number of electrons per unit 
cell, should be integrated over a surface, on which these 
N bands are separated by an energy gap from the higher 
energy ones. Since the crossing of the A'th and N + 1th 
bands occurs at a point (Weyl point), such a surface en¬ 
closing the type-II Weyl point always exists. The calcula¬ 
tion becomes equivalent to the usual integration over the 
occupied bands, if one allows for a /c-dependent chemi¬ 
cal potential, located in between the A'-th and N -f 1-th 
bands at each /c-point on the surface. 

Following this prescription and using a Wannier 
function-based tight-binding model [49, 50] described be¬ 
low, Bloch states were calculated on the spheres around 
the gapless points. The Hamiltonian remained gapped in 
the above sense (A = 72) everywhere on these spheres, 
schematically illustrated as ii ^2 in Fig- 3(a) of the main 
text. The corresponding flux of Berry curvature was com¬ 
puted by discretizing a closed sphere, parametrized by 
angles 0 and 0, into ID-loops, as shown in Fig. 14(a). 
When traversed around a loop, the occupied Bloch states 
accumulate a Berry phase, the trace of which is com¬ 
puted for each of these loops Oi using the methods of 
Refs. [21, 22, 51]. 

These Berry phases correspond to the average position 
of charge [52], associated with the bands below the gap, 
on the loop, {(p){0i). Since ID loops cover a closed sur¬ 
face, the center of charge {(p) can only shift by an integer 
multiple of 27r when 0 varies from 0 to tt. This multiple 
is equal to the monopole charge enclosed and gives the 
chirality of the WP enclosed. Similar considerations are 
used when computing Chern numbers of insulators, with 
the only difference that the BZ in that case is a torus, 
rather than a sphere. The result of such a calculation is 
equivalent to taking the surface integral of the Berry cur- 
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Figure 15. Schematic illustration of the effect of uniaxial com¬ 
pressive strain along the ^-direction. Left panel: at the strain 
of ~2% four Weyl points meet in pairs on the mirror plane 
and annihilate. Right panel: only four Weyl points are left in 
the kz — ^ plane. 

vature over a closed surface. A more detailed account of 
this type of calculations, along with the rigorous deriva¬ 
tion can be found in Ref. [? ]. 

The results obtained for WTe 2 is shown in Fig. 14(b) 
for the sphere of Fig. 3(a) of the main text. The charge 
center shifts downwards, corresponding to Chern number 
C = — 1, thus proving the existance of a WP of negative 
chirality inside li. For the sphere £ 2 ^ enclosing the second 
Weyl point, the Chern number is C = +1, and hence the 
chirality of this point is positive. 


Appendix G: Strain effects in WTe 2 

As mentioned in the main text, the hydrostatic com¬ 
pressive strain (applied pressure) of 2% increases the sep¬ 
aration of the Weyl nodes from 0.7% to 4% of the re¬ 
ciprocal lattice constant separation of the nearest Weyl 
points. Other strains were also studied. Here we present 
the summary of the results obtained. 

Compressive uniaxial strain along the x-direction also 
increases the separation between the nearest Weyl points. 
At 1% strain the separation is 2.2% of |G 2 |- Further in¬ 
crease of this strain makes one of the two points hit the 
mirror plane (at ~2% strain), where it is annihilated with 
its mirror image of opposite chirality, and thus only four 
Weyl points are left in the kz = ^ plane, all of them being 
of type-IL This situation, illustrated if Fig. 15, possibly 
realizes the simplest situation in time-reversal symmet¬ 
ric Weyl semimetals with the minimal number of Weyl 
points. 

The situation changes for the uniaxial compressive 
strain in the ^-direction. While the the neighboring 
Weyls still move further from each other (4.3% of IG 2 I 
for a 2% strain), they both move further away from the 
mirror plane, so the scenario of only four Weyl points 
left is not realized in this case. However, we find that 
this strain drives a phase transition from type-H to type- 
I Weyl points for four out of eight points. Those, that 
are furthest from the mirror plane become type-I at ~2% 
strain as illustrated in Fig. 16. 

Stretching the crystal in both x and ^ directions leads 



Figure 16. Illustration of type-H to type-I Weyl point transi¬ 
tion under the compressive uniaxial strain in the ^-direction. 
Left panel: no strain. All points are of type-H. Right panel: 
2% strain. The points at higher energy were checked in 3D k- 
space to be of type-I. The points move further apart in energy, 
one of them appearing very close to Ey. 

to the mutual annihilation of the neighboring Weyl points 
and drives the transition from a topological semimetal 
with eight Weyl points to the trivial semimetal with no 
Weyl points at all. For the uniaxial ^-strain the situa¬ 
tion is reversed. Stretching moves the neighboring Weyl 
points away from each other (2.5% IG 2 I for a 2% stretch) 
and for reasonable strains the number of Weyl points does 
not change, while compression in this direction leads to 
annihilation of all points. 

We also considered a strain in the [lll]-direction, 
which breaks all the symmetries but time-reversal. The 
Weyl points survive at small strains, but they move off 
the high-symmetry kz = ^ plane. This should be ex¬ 
pected, since the rotational G 2 symmetry, which in com¬ 
bination with time-reversal allowed the Weyl points to 
appear on the high-symmetry plane, is now broken, and 
a stable Weyl point should appear at arbitrary /c-points. 

Appendix H: Tight-binding models 

To compute chiralities (Chern numbers) of Weyl 
points, as well as to interpolate between no SOC and 
full SOC states, tight-binding models were derived with 
and without SOC. These tight-binding models were ob¬ 
tained using Wannier interpolation [49, 50]. Bloch states 
were projected onto all W d-states and all Te p-states. 
We used a specially symmetrized model without SOC 
that has the same Weyl points distribution as obtained 
from first-principles, and interpolated it to the one de¬ 
rived from a fully relativistic calculation. 

To get the chiralities of Weyl points for the full SOC 
case, we used an 88-band tight binding model derived 
from the full SOC first-principles calculation. 


Appendix I: Surface states 

Here we provide an illustration to support our asser¬ 
tion, made in the main text, about the possibility to move 
the Fermi arc that crosses the hole pocket into the con¬ 
tinuum of the bulk states. Fig. 17 shows a larger scale 
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Figure 17. Spectral function illustrates the surface Fermi sur¬ 
face. The surface state, indicated with arrows, emerges in the 
electron pocket (shown with the black contour), crosses the 
hole pocket (shown with the white contour), and merges back 
into the electron pocket (by mirror symmetry). 


illustration, where it can be seen that this surface state 
emerges from the electron pocket and merges back into 
it. Note that unlike the illustration of Fig. 4(b) of the 
main text, the small electron pocket becomes part of the 
larger electron pocket in this illustration. This is due to 
small deviations lOmeV) of our tight-binding model 
band structure from its ab initio counterpart used for 
illustrations in Fig. 4 of the main text. 


Appendix J: Additional Remarks 


Matplotlib [53], Mayavi [54] and VESTA [55] software 
packages were used to create some of the illustrations. 
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